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Abstract. In this paper we construct a global, continuous flow of solutions to the Camassa-Holm 
equation on the entire space H 1 . Our solutions are conservative, in the sense that the total energy 
J (u 2 + u 2 ) dx remains a.e. constant in time. Our new approach is based on a distance functional 
J(u,v), defined in terms of an optimal transportation problem, which satisfies -^J(u(t),v(t)) < 
k ■ J(u(t),v(t)) for every couple of solutions. Using this new distance functional, we can construct 
arbitrary solutions as the uniform limit of multi-peakon solutions, and prove a general uniqueness 
result. 



1 - Introduction 

The Camassa-Holm equation can be written as a scalar conservation law with an additional 
integro-differential term: 

iH + {u 2 /2) x + P x = 0, (1.1) 



where P is defined as a convolution: 



P = L-W* ( u 2 + 4- )• (1-2) 
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For the physical motivations of this equation we refer to [CH], [CM1], [CM2], [J]. Earlier results 
on the existence and uniqueness of solutions can be found in [XZ1], [XZ2]. One can regard (1.1) 
as an evolution equation on a space of absolutely continuous functions with derivatives u x £ L 2 . 
In the smooth case, differentiating (1.1) w.r.t. x one obtains 

.2 f..2 u2 



u xt + uu xx + u x - [u +Yj + P = - (I- 3 ) 

Multiplying (1.1) by u and (1.3) by u x we obtain the two conservation laws with source term 

u 2 \ ( u 3 
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As a consequence, for regular solutions the total energy 

E(t) = J [u 2 (t,x) +u 2 x (t,x)]dx (1.6) 

remains constant in time. 

As in the case of scalar conservation laws, by the strong nonlinearity of the equations, solutions 
with smooth initial data can lose regularity in finite time. For the Camassa-Holm equation (1.1), 
however, the uniform bound on \\u x ||l 2 guarantees that only the L°° norm of the gradient can blow 
up, while the solution u itself remains Holder continuous at all times. 

In order to construct global in time solutions, two main approaches have recently been intro- 
duced. On one hand, one can add a small diffusion term in the right hand side of (1.1), and recover 
solutions of the original equations as a vanishing viscosity limit [CHK1, CHK2]. An alternative 
technique, developed in [BC2], relies on a new set of independent and dependent variables, specif- 
ically designed with the aim of "resolving" all singularities. In terms of these new variables, the 
solution to the Cauchy problem becomes regular for all times, and can be obtained as the unique 
fixed point of a contractive transformation. 

In the present paper, we implement yet another approach to the Camassa-Holm equation. As 
a starting point we consider all multi-peakon solutions, of the form 

N 

u(t,x) = J2P^ e ~ lX ~ qiit)l - ( L7 ) 

i=l 

These are obtained by solving the system of O.D.E's 

1 j (1.8) 

p i = ^2p i p j sign(q i -q j )e l?l q ^ . 

It is well known that this can be written in hamiltonian form: 

8 , . 
Qi = -K-H(p,q), 

dqi 



If all the coefficients Pi are initially positive, then they remain positive and bounded for all times. 
The solution u = u(t, x) is thus uniformly Lipschitz continuous. We stress, however, that here we 
are not making any assumption about the signs of the pi. In a typical situation, two peakons can 
cross at a finite time r. As t — > r— their strengths Pi,Pj and positions qi, qj will satisfy 

Pi{t) ^+oo, pj(t) -y-oo, pi(t)+pj(t) ->p, (1.9) 

q l {t)^q, qj(t)^q, qi (t) < qj(t) for t < r , (1.10) 
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for some p,q 6 JR. Moreover, ||^x(i)|| Loo - > °°- in this case, we will show that there exists a unique 
way to extend the multi-peakon solution beyond the interaction time, so that the total energy is 
conserved. 

Having constructed a set of "multi-peakon solutions", our main goal is to show that these 
solutions form a continuous semigroup, whose domain is dense in the space H 1 ^!!!). Taking the 
unique continuous extension, we thus obtain a continuous semigroup of solutions of (1.1), defined 
on the entire space H 1 . 

The particular metric, used to derive the Lipschitz continuity of the semigroup, is the most 
novel feature of our paper. Indeed, one easily checks that the flow map &t '■ u(0) >— ► u(t) cannot 
be continuous as a map from H 1 into itself, or from L 2 into itself. Distances defined in terms of 
convex norms perform well in connection with linear problems, but occasionally fail when nonlinear 
features become dominant. In the present setting, we construct a new distance J(u,v) between 
functions u, v £ H , defined by a problem of optimal transportation. Roughly speaking, J(u,v) 
will be the minimum cost in order to transport the mass distribution with density 1 + u 2 located 
on the graph of u onto the mass distribution with density 1 + v 2 located on the graph of v. See 
Section 3 for details. With this definition of distance, our main result shows that 

< C-J(u(t), v(t)) (1.11) 

for some constant C and any couple of multi-peakon solutions u, v. Moreover, J(u n ,u) — > implies 
the uniform — * 0. The distance functional J thus provides the ideal tool 
to measure continuous dependence on the initial data for solutions to the Camassa-Holm equation. 
Earlier applications of distances defined in terms of optimal transportation problems can be found 
in the monograph [V]. 

The issue of uniqueness of solutions must here be discussed in greater detail. For a multi- 
peakon solution, as long as all coefficients pi remain bounded, the solution to the system of 
ODE's (1.8) is clearly unique. For each time t, call [it the measure having density u 2 {i) + u 2 ,^) 
w.r.t. Lebesgue measure. Consider a time r where a positive and a negative peakon collide, accord- 
ing to (1.9)-(1.10). As t — ► t— , we have the weak convergence jit — 1 [i T for some positive measure 
fi T which typically contains a Dirac mass at the point q. By energy conservation, we thus have 

J [u 2 (t, x) + u 2 (t,x)] dx + fi T ({q}) = lim J [u 2 (t, x) + u 2 , (r, x)] dx = E(r—) . 

There are now two natural ways to prolong the multi-peakon solution beyond time r: a conservative 
solution, such that 

E(t) = J [u 2 (t,x) + u 2 x (t,x)]dx = E{t-) t>T, 

or a dissipative solution, where all the energy concentrated at the point q is lost. In this case 

E{t) = j [u 2 (t,x) + u 2 x (t,x)]dx = E{t-)~ fi T ({q}) t>T. 

For t > t, the dissipative solution is obtained by simply replacing the two peakons Pi,Pj with one 
single peakon of strength p, located at x = q . On the other hand, as we will show in Section 2, the 
conservative solution contains two peakons emerging from the point q. As t — » r+, their strengths 
and positions satisfy again (1.9), while (1.10) is replaced by 

fc(t)->g, Qj(t)^q, Qi{t)>qj{t) for t>r. (1.12) 



-J{u(t), v(t)) 
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The vanishing viscosity approach in [CHK1, CHK2] singles out the dissipative solutions. These 
can also be characterized by the Oleinik type estimate 

u x (t,x) < C(l + t" 1 ), 

valid for t > at a.e. x 6 JR. On the other hand, the coordinate transformation approach in [BC2] 
and the present one, based on optimal transport metrics, appear to be well suited for the study of 
both conservative and dissipative solutions. 

In the present paper we focus on conservative solutions to the C-H equation in the spatially 
periodic case. The assumption of spatial periodicity allows us to concentrate on the heart of the 
matter, i.e. the uniqueness and stability of solutions beyond the time of singularity formation. 
It will spare us some technicalities, such as the analysis of the tail decay of u, u x as x — » ±00. 
The construction of an optimal transportation metric in connection with dissipative solutions to 
the C-H equation is significantly different, and will be carried out in a forthcoming paper. The 
main ingredients can already be found in the paper [BC1], devoted to dissipative solutions of the 
Hunter-Saxton equation. 

As initial data, we take 

u(0,x) = u(x), (1.13) 

with u in the space Hp ei of of periodic, absolutely continuous functions u with derivative u x € Lp er . 
To fix the ideas, we assume that the period is 1, so that 



u(x + 1) = u(x) x <G M . 

On iT* we shall use the norm 



,1 2 ,1 2 x 1/2 

|it||zfi er = I / dx + \u x (x)\ dx 



Our main results can be stated as follows. 



Theorem 1. For each initial data u 6 H^ ex , there exists a solution u(-) of the Cauchy problem 
(1.1), (1.13). Namely, the map 1 1— > u(t) is Lipschitz continuous from M into Lp er , satisfies (1.13) 
at time t = 0, and the identity 

—u = -uu x - P x (1.14) 

is satisfied as an equality between elements in Lp er at a.e. time t G M. This same map t <— ► u(t) 
is continuously differentiable from M into L^ er and satisfies (1.14) at a - e - time t G JR, for all 
p £ [1,2[. The above solution is conservative in the sense that, for a.e. t G M, 

E(t)= [ [u 2 (t,x)+u 2 x (t,x)]dx = E u = [ [u 2 (x)+u 2 x (x)]dx. (1.15) 
Jo Jo 



Theorem 2. Conservative solutions to (1.1) can be constructed so that they constitute a continuous 
flow Namely, there exists a distance functional J on H^ er such that 

q ■ \\u-v\\ Her < J(u,v) <C- ||u-u||jri er (1.16) 
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for all u, v £ H^ er and some constant C uniformly valid on bounded sets of H^ ei . Moreover, for 
any two solutions u(t) = $ t u, v(t) = $ t v of (1.1), the map t i— >• j(u(t), v(t)) satisfies 

J(u(t), u) < d ■ \t\ , (1.17) 

J(u(t),v(t)) < J(u,v)-e C2W (1.18) 
for a.e. t € M and constants C\,C<i, uniformly valid as u,v range on bounded sets of H^ er . 

Somewhat surprisingly, all the properties stated in Theorem 1 are still not strong enough to 
single out a unique solution. To achieve uniqueness, an additional condition is needed. 

Theorem 3. Conservative solutions t ^ u(t) of (1.1) can be constructed with the following 
additional property: 

For each t € M, call /i t the absolutely continuous measure having density u 2 + u x w.r.t. Lebesgue 
measure. Then, by possibly redefining fj, t on a set of times of measure zero, the map t i— > fj, t 
is continuous w.r.t. the topology of weak convergence of measures. It provides a measure-valued 
solution to the conservation law 

w t + (uw) x = (u 3 - 2uP u ) x . (1.19) 



The solution of the Cauchy problem (1.1), (1-13) satisfying the properties stated in Theorem 
1 and this additional condition is unique. 

The plan of the paper is as follows. In Section 2 we derive some elementary properties of 
multi-peakon solutions and show that any initial data can be approximated in Hp er by a finite sum 
of peakons. In Section 3 we introduce our distance functional J(u, v) and study its relations with 
other distances defined by Sobolev norms. The continuity of the flow (1.1), together with the key 
estimates (1.17)-(1.18) are then proved in the following two sections. The proofs of Theorems 1 
and 2 are completed in Section 6. The uniqueness result stated in Theorem 3 is proved in Section 7. 
As a corollary, we also show that in a multi-peakon solution the only possible interactions involve 
exactly two peakons: one positive and one negative. In particular, no triple interactions can ever 
occur. 



2 - Multipeakon solutions 

By a periodic peakon we mean a function of the form 

u(x)=p X (x-q), X {x) = E e ~ lX ~ nl • 

Observe that the periodic function \ satisfies 

x{-%) = x(x) = x(x + 1) xeR, 

i ^ e^ + e 1 "* [n nl (2-2) 

X{x) = — x G [0, 1] . 

e — 1 



We begin this section by observing that any periodic initial data can be approximated by a 
periodic multi-peakon. 

Lemma 1. Let f G H^ er . Then for any e > there exists periodic multi-peakon g, of the form 

N N 

9{x) = J> e- lx - q *- nl = J> X{x ~ ft) (2.3) 

i=l n£Z i=l 

stic/i that 

ll/-<?ll^ er <£• 

Proof. By taking a suitable mollification, we can approximate / with a periodic function / 6 C°°, 
so that 

||/-/|| H i er <e/2. (2.4) 

Next, we observe that 




where So denotes the Dirac distribution concentrating a unit mass at the origin. We can thus write 
/ as a convolution: 

/=*'/=K e " H -^ e " w )* /=e " H *( £ ^)- 

The above integral can now be approximated with a Riemann sum 

< \ i \ r /N M - im j 

9{x) = ]>PiX{x ~ 1i) , Pi= o dy. 

i=l J{i-1)/N 2 

Choosing iV sufficiently large we obtain ||/ — g\\H^ er < e/2. Together with (2.4) this yields the 
result. □ 

Next, we show how to construct a unique conservative solution, for multi-peakon initial data. 
As long as the locations ft of the peakons remain distinct, this can be obtained by solving the 
Hamiltonian system of O.D.E's (1.8). However, at a time r where two or more peakons interact, the 
corresponding strengths pi become unbounded. A suitable transformation of variables is needed, 
in order to resolve the singularity and uniquely extend the solution beyond the interaction time. 

Lemma 2. Let u be any periodic, multi-peakon initial data. Then the Cauchy problem (1.1), 
(1.13) has a global, conservative multi-peakon solution defined for all t G JR. The set T of times 
where two or more peakons interact is at most countable. Moreover, for all t G" 1, the energy 
conservation (1.15) holds. 
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Proof. The solution can be uniquely constructed by solving the hamiltonian system (1.8), up to 
the first time r where two or more peakons interact. We now show that there exists a unique way 
to prolong the solution for t > r, in terms of two outgoing peakons. To fix the ideas, call 



q= lim qi(t) i = l,...,k, 

the place where the interaction occurs, and let pi(t), . . . ,Pk(t) be the strengths of the interacting 
peakons. Later in this paper we will show that only the case k = 2 can actually occur, but at this 
stage we need to consider the more general case. We observe that the strengths Pk+i, ■■■ ,Pn of 
the peakons not involved in the interaction remain continuous at time r. Moreover, by (1.8) there 
exists the limit 

k 

p= lim Vpi(i). 

t^T— *■ » 

i=l 

We can thus write 



N N 

\mi_J2Pi(t)e~ lx ~ qi{t)l = pe~\ x -^+ ^ Pi (r) e"!*-*™ 

i=l i=k + \ 



For t > r, we shall prolong the solution with two peakons emerging from the point q. The strength 
of these two peakons will be uniquely determined by the requirement of energy conservation (1.15). 

Call £ + (i) respectively the position of the smallest and largest characteristic curves 

passing through the point (t, q), namely 



£-(*) = min {£(*); £(r)=g, £(s) = u(s,£(s)) for all s G [r - h,r + h]} , 
£ + (t) = maxj^(t) ; £(r) = g, f(s) = «(a,^(s)) for all s G [t - h,r + h]} . 



(2.5) 



Moreover, define 

/■€+(*) 

e (r,g) = t nm / U x (t,x)dx. (2-6) 



The existence of this limit follows from the balance law (1.5). This describes how much energy is 
concentrated at the interaction point. 

For t > t the solution will contain the peakons Pk+i, ■ ■ ■ ,Pn, located at qk+i, ■ ■ ■ ,qN, together 
with the two outgoing peakons P i,P2, located at q\ < q2- The behavior of P i,qi for i G {k + 
1, . . . , N} is still described by a system of O.D.E's as in (1.8). However, to describe the evolution 
of PiiP2,Qi,Q2 one has to use a different set of variables, resolving the singularity occurring at 
(r, q). As t — > r+ we expect (1.9), (1.12) to hold. To devise a suitable set of rescaled variables, we 
observe that, by (1.3), 

±u x (t,m) = -\ul(t,at)) + [n 2 ~ P] (2.7) 

along any characteristic curve t i-> £(t) emerging from the point q. Since u,P remain uniformly 
bounded, one has 

2 

u x (t,x)&- t>T, x G [q^t) , q 2 (t)] . 
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The total amount of energy concentrated in the interval between the two peakons is given by 



/ [u 2 (t, x) + u x (t, x)] dx 

Jqi(t) 



u(q 2 )-u(q 1 )\ . . 

-(92-91) 

92-91 / 
(Pi - Pi) 2 (92 - 9i ) ~ e (T)? ) . 



(P2 -Pi)(l - e-l^-^l) 



92 - 9i 



The previous heuristic analysis suggests that, in order to resolve the singularities, we should work 
with the variables 

z=p!+p 2 , w = 2arctan(>2 ~ Pi) , V = q 2 + qi , ( = (p 2 ~ Pi f (q 2 ~ qi) , 

together with pk+i , • • • , Pn , 9fc+i , • • • , 9jv- To simplify the following calculations we here 
assume < q\ < q 2 < qk+i < ••• < 9jv < 1 5 which is not restrictive. 

Let x defined in (2.1) and x{ x ) = ~ e "t^ — , x € [0, 1]. From the original system of equations 
(1.8) it follows 



z = cosh 



2tan 2 f 



x N oinh ( i ^ N 

\ 2 w ( V\ llm l2tan 2 f ) ^ / 77\ 

\ zcos ^ PjX{Qj -gj IXI PiX^-jJ 

' 7=fe+l " tan f 7 = fc+l 



2"/ 



c 



+ 2 cosh 



7 j=k+l 
\ 

Z 



2tan 2 ^, 

2 / j=fc+i 



J] Pj-x(?,--|) 



+ 2 sinh 



2 tan 2 f 



sinu; 



N 

Y p i*{ 



7] = Z 



x(o) + x 



tan 2 § 



j=fc+i 



+ cosh 



*"2 



2tan z 



-) 2 



■ _X(0)-X(t5^ 



1 

tan 2 ¥ 



^^ X{ tan 2 H J tan 1 



sinh ( 771—^ 

V 2 tan 



2 w 
2 



+ 2- 



tanf 



cosh 



2 / 2 

AT 



1 

tan 2 4^ 



N 

Y Pjx(ij-l) 



j=k+l 



Z 



2 tan 2 f / tan „ 



N , . s JV 

- ;— t 1 * 2 / j = k-\-i 



Pi =Pi 



cosh 



N 



2 tan 2 f 



^X ( Qi ~ |) + sinh 



tan— x ( qi - - J 



2i an 2 f 7 2 



+ Pi Y p J si S n (9i - 9j)x (l9i - 9j I) 



j=k+i 



qi = cosh 



2 tan 2 f 



( ft " I) 



+ sinh 



2 tan' 



\ N 

) tan (ft ~^)+ Y p i x (I* ~ q i 

2 / " j=k+l 



with initial data 



z(r)=p, w {t)=tt, T](r) = 2q, ((t) = e {r4) , 
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Pi(r) = limpet), 9i(t) = lim q^t) i = k + l,...,N. 

For the above system of O.D.E's, a direct inspection reveals that the right hand side can be 
extended by continuity also at the value w = tt, because all singularities are removable. This 
continuous extension is actually smooth, in a neighborhood of the initial data. Therefore, our 
Cauchy problem has a unique local solution. This provides a multi-peakon solution defined on 
some interval of the form [r, t'[ , up to the next interaction time. 

The case where two or more groups of peakons interact exactly at the same time r, but at 
different locations within the interval [0,1], can be treated in exactly the same way. Since the 
total number of peakons (on a unit interval in the x-variable) does not increase, it is clear that the 
number of interaction times is at most countable. The solution can thus be extended to all times 
t > 0, conserving its total energy. □ 



3 - A distance functional 

In this section we shall construct a functional J(u, v) which controls the distance between two 
solutions of the equation (1.1). All functions and measures on M are assumed to be periodic with 
period 1. Let T be the unit circle, so that T = [0, 2tt] with the endpoints and 2tt identified. The 
distance \9 — 9'\* between two points 9, 9' G T is defined as the smaller between the lengths of the 
two arcs connecting 9 with 9' (one clockwise, the other counterclockwise). We now consider the 
product space 

X = M x M x T 

with distance 

(^(x, u, w), (x,u,w)^ = (^\x — x\ + \u — u\ + | io — w\*^J A 1 , (3-1) 

where a A b = min{a, b}. Let M(X) be the space of all Radon measures on X which are 1-periodic 
w.r.t. the x- variable. To each 1-periodic function u G H^ er we now associate the positive measure 
a u G M{X) defined as 

° U {A)=L (l + ul(x))dx (3.2) 

J \x£M: (x, u(x), 2 arctan u x (x) )dAf 

for every Borel set A C M 2 x T . Notice that the total mass of a^ u '^ over one period is 

tj u ([0, 1] x M x T) = 1 + f u 2 x (x) dx . 

Jo 

On this family of positive, 1-periodic Radon measures, we now introduce a kind of Kantorovich 
distance, related to an optimal transportation problem. Given the two measures a u and a u , their 
distance J(u, u) is defined as follows. 

Call T the family of all strictly increasing absolutely continuous maps tj; : 1R M which have 
an absolutely continuous inverse and satisfy 

ip(x + n) = n + tp(x) for every n G /Z . (3.3) 
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(3.4) 



For a given -0 € T, we define the 1-periodic, measurable functions 0i, 02 : M i-> [0, 1] by setting 
0!(x) = sup jfle [0,1]; 0-(l + ^(x)) < (l + ^(^(x)))0'(x 
2 (x) =sup jfl € [0,1]; 1 + «2(ar) > 0- (l + u\ (ip(x))^J ip' (a 

Observe that the above definitions imply max. {<j>i(x) , 2 {%)} = 1 together with 

^(l + u^x)) =0 2 (0(x)) (l + u 2 x {^(x)))^(x) (3.5) 



for a.e. x G JR. We now define 



J^(u,u) = y ^(x, -u(x), 2 arctanu^x)) , (ip(x), u(ip(x)), 2 arctan-u ;E (-0(x))^ 

■0i(x) (1 + ^0*0)^ (3.6) 



(l + n 2 (x)) _ (l + S 2 Wa;))) ^ (x) 



<ix . 



Of course, the integral is always computed over one period. Observe that x *— > ^A(x) can be regarded 
as a transportation plan, in order to transport the measure a u onto the measure a u . Since these 
two positive measures need not have the same total mass, we allow the presence of some excess 
mass, not transferred from one place to the other. The penalty for this excess mass is given by 
the second integral in (3.6). The factor 0i < 1 in the first integral indicates the percentage of the 
mass which is actually transported. Integrating (3.5) over one period, we find 

/ 1 (x)(l + ^(x))dx = [ (f) 2 {y){l + ul{y))dy. 
■Jo Jo 

We can thus transport the measure 0i a u onto 02 cr u by a map Vl/ : (x, u(x) arctanu^x)) i— ► 
(y, u(y), arctann^y)), with y = tp(x). The associated cost is given by the first integral in (3.6). 
Notice that in this case the measure 02 a u is obtained as the push- forward of the measure 0i a u . We 
recall that the push-forward of a measure a by a mapping \& is defined as (ty#a)(A) = a(^~ 1 (A)) 
for every measurable set A. Here ^~ 1 (A) = {z ; *£(z) € A}. 

Our distance functional J is now obtained by optimizing over all transportation plans, namely 

J(u,u) = inf J^(u,u). (3.7) 

To check that (3.7) actually defines a distance, let u,v,w 6 H 1 (M) be given. 

1. Choosing ip(x) = x, so that 0i(x) = 2 (:c) = 1 5 we immediately see that J(u,u) = 0. Moreover, 
if J(u, u) = 0, then by the definition of d^ we have u = u. 

2. Given tp <G T, define tp = -0 -1 , so that 0i = 02, 02 = 0i- This yields 
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Hence J(u,u) = J(u,u). 

3. Finally, to prove the triangle inequality, let ip* , : 2R i — >• 1R be two increasing diffeomorphisms 
satisfying (3.3), and let <p\ , <f)\ , <f)\ , <t>\ '■ ^ l— * [0, 1] be the corresponding functions, defined as in 
(3.4). We now consider the composition -0 = o -0 b and define the functions 4>i, (f>2 according to 
(3.4). Observing that 

MX) >0i(x)-^(V b (x)), 

MM*)) = h(¥{M(x))) > M^( x )) ■ 4(¥W(x))) , 

and recalling that the distance at (3.1) is always < 1, we conclude 

(u, v ) + (v, w) . 

This implies the triangle inequality J(u,v) + J(v,w) > J(u,w). □ 
In the remainder of this section we study the relations between our distance functional J and 
the distances determined by various norms. 



Lemma 3. For any u, v G H^ er one has 

u-v\\ L i < J(u,v) < C ■ \\u-v\\ H i , (3.8) 



'per 

1 

c 

with a constant C uniformly valid on bounded subsets of H^ er . 

Proof. We shall use the elementary bound 

| arctana — arctanfrj • a 2 < 47r(|a| + |6|) \a — b\ , (3.9) 
valid for all a,b G JR. In connection with the identity mapping tp(x) = x we now compute 

J^(u,v) < J ^\u(x) — v(x)\ + 2| arctanu x — arctant^lj (1 + u 2 ) dx + J \u 2 x — v 2 \dx 

< \\ u — ||1 + ^IIl 1 + (87r + 1) / \u x + v x \ \u x — v x \ dx 

Jo 

< (8tt + 3) (l + \\u\\ H i + |M|#i) • \\u - v\\ H i , 
proving the second inequality in (3.8). 
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To achieve the first inequality, choose any ip G T. For x G [0, 1], call j x the segment joining 
the point P x = (x,u(x)^ with Q x = (Y>(;c), v(tp(x))^j . Clearly, the union of all these segments 
covers the region between the graphs of u and v. Moving the base point from x to x + dx, the 
corresponding segments sweep an infinitesimal area dA estimated by (fig. 1) 



\dA\ < \P X -Q x \ ■ (\dP x \ + \dQ x \) 

1/2 

< (|x-V(x)| 2 + |n(x)-t;(^(x))| 2 + ) • [{l + u 2 x ) 1 / 2 dx + {l + vl) 1 ' 2 ^'{x)dx 
Integrating over one period we obtain 

u{x) — v(x) | dx 

< / (\ x ~ *(t)(x)\ + \u(x) - t# (t) (*))|) • [(1 + ul(x)) + (1 + vl^{x))) 1,2 ^{ t) {x) 

< (2+ \\u\\ H i + \\v\\ H i) ■ \j^(u,v) + J^~\v,u)\ < C-J(u,v), 

completing the proof of (3.8). 



dx 



(3.10) 
□ 



Lemma 4. Let (« n ) n >i be a Cauchy sequence for the distance J, uniformly bounded in the Hp er 
norm. Then 

(i) There exists a limit function u € H^ ex such that u n — > u in L°° and the sequence of derivatives 
u n ,x converges to u x in L^ er , for 1 < p < 2. 

(ii) Let [i n be the absolutely continuous measure having density u\ x w.r.t. Lebesgue measure. Then 
one has the weak convergence fi n — ^ \i, for some measure fx whose absolutely continuous part 
has density u\. 



Proof. 1. By Lemma 3 we already know the convergence u n — > u, for some limit function 
u G Lp er . By a Sobolev embedding theorem, all functions u n ,u are uniformly Holder continuous. 
This implies \\u n — - ► 0. 

To establish the convergence of derivatives, we first show that the sequence of functions 

v n = exp{2iarctanti nja ;} 

is compact in Lp er . Indeed, fix e > 0. Then there exists N such that J(u m , u n ) < e for m,n > N. 
We can now approximate un in H^ ev with a piecewise affine function un such that J(un,Un) < e. 
By assumption, choosing suitable transport maps ip n we obtain 

dx < 2J(u n ,UN) < 4e 

for all n> N. We now observe that all functions x i— > exp \2i arctanu7\r,a;(V'n(a ; ))} are uniformly 
bounded, piecewise constant with the same number of jumps: namely, the number of subintervals 
on which un is affine. The set of all such functions is compact in L per . This argument shows 
that the sequence v n = exp{2i arctan-u nja ;} eventually remains in an e-neighborhood of a compact 



exp { 2i arctan u n ^ x (x) } — exp { 2i arctan u n, x (VVi ( x ) ) } 
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subset of L per . Since e > can be taken arbitrarily small, by possibly choosing a subsequence we 
obtain the strong convergence v n — > v for some v G L per . 

2. From the uniform H 1 bounds and the L 1 convergence of the functions v n , we now derive the 
L p convergence of the derivatives. For a given e > 0, define 

M = sup \\u n \\m ei , A n = \x G [0, 1] ; \u n>x (x)\ > M/e\ . 

The above definitions imply 

We now have 

||^rrc,x 'Wn^llLP 5: 



meas(A„) < e 2 



(3-11) 



,l]\(A„UA m ) 



= /l + I 2 ■ 



n (2-p)/2p r 



1 • dx 



< £ (2-p)/p . 2 Af . 

Next, choosing a constant C £ such that 



/ (|«m,x| + \u n ,x\) dx 

J A m. U -A n. 



(3.12) 
(3.13) 



we obtain 



a 2i arctan a arctan b 



h<c £ 



> C e \a-b\ 



D 2i arctan u miX _ arctan u„ lX \P ^ 



whenever \a\, \b\ < M/e , 



(3.14) 



Taking e > small, we can make the right hand side of (3.13) as small as we like. On the other 
hand, choosing a subsequence such that v v = e 2iarctanu ^ converges in L per , the right hand side 
of (3.154) approaches zero. Hence, for this subsequence, 

limsup \\u m<x -u„ ia; || L P er =0. 

m,n— roo 

Since u n — ► u uniformly, in this case we must have 



\U>r).. X ^ 



0. 



(3.15) 



We now observe that from any subsequence we can extract a further subsequence for which (3.15) 
holds. Therefore, the whole sequence (« n ,x) n >i converges to u x in L^ er . 

3. To establish (ii), we consider the sequence of measures having density 1 + u 2 x w.r.t. Lebesgue 
measure. This sequence converges weakly, because our distance functional is stronger than the 
Kantorovich-Waserstein metric which induces the topology of weak convergence on spaces of mea- 
sures. Therefore, \i n — 1 \i for some positive measure \i. 

Since the sequence 1 + u n ^ x converges to l + u x in L per , by possibly choosing a subsequence we 
achieve the pointwise convergence u nyX {x) — > u x (x), for a.e. x G [0, 1]. For any e > 0, by Egorov's 
theorem we have the uniform convergence u n>x (x) — > u x (x) for all x G [0, 1] \ V s , for some set with 
meas(T4) < £■ Since e > can be taken arbitrarily small, this shows that the absolutely continuous 
part of the measure ji has density v? + v? x w.r.t. Lebesgue measure. □ 
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4 - Continuity in time of the distance functional 

Here and in the next section we examine how the distance functional «/(•,•) evolves in time, 
in connection with multi-peakon solutions of the Camassa-Holm equation (1.1). We first provide 
estimates valid on a time interval where no peakon interactions occur. Then we show that the 
distance functional is continuous across times of interaction. Since the number of peakons is 
locally finite, this will suffice to derive the basic estimates (1.1 7)- (1.18), in the case of multi-peakon 
solutions. 

Lemma 5. Let t u(t) G H^ er be a multi-peakon solution of (1.1). Assume that no peakon 
interactions occur within the interval [0, r] . Then 

J(u(s),u(s')) <C -\s- s'\, for alls, s' g[0,t], (4.1) 

for some constant C, uniformly valid as u ranges on bounded subsets of H^ ex . 

Proof. Assume < s < s' < r. By the assumptions, the solution u = u(t, x) remains uniformly 
Lipschitz continuous on the time interval [0,r]. Therefore, for each s G [0,r] and x G M, the 
Cauchy problem 

±C(t)=u(t,m), C(s)=x, (4.2) 

determines a unique characteristic curve t i— > £(t;s,x) passing through the point (s,x). Given 
s' 6 [0, r], we can thus define a transportation plan by setting 

i>(x) = as';s,x). (4.3) 

Of course, moving mass along the characteristics is the most natural thing to do. We then choose 
(f>i,4>2 to be as large as possible, according to (3.4). Namely: 

^i(x) = sup jfl G [0, 1] ; 9 ■ (l + u 2 x (s, xfj < i//(x) • (l + u 2 x (s' , y>(x))) j , 

<h{x) = SU P {^e [0,1]; 6-(l + ul(s', #c))) tf(x) < l + i&s.x)}. 
The cost of this plan is bounded by 

J^(u(s), u(s'))<J {\x-Z(s';s,x)\ + \u(s,x)-u(s',Z(s';s,x))\ 

+ |2 arctan Ua;(s, x) — 2 arctanw x (s', £(s'; s, x))|^| (l + u x (s, x)) dx 

+ I {l + ul{s,x))dx 

Jo 

+ J (1 - MM*))) (l + V^))) V'(ar) . 
To estimate the right hand side of (4.4), we first observe that, for all u G H^ eT , 

< / |ii(x)|(fx+ / |« x (ar)|da; 
Jo Jo 

< ||u|| L 2 + ||-u x || L 2 < 2\\u\\ H i = 2(E U ) 1/2 . 

ii ii ii ii 1 1 1 1 per \ / 
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(4.4) 



(4.5) 



Using (4.5) in (4.2) we obtain 

\as)-as')\<2m i / 2 .\ s - s '\. 

Next, from the definition of the source term P at (1.2) it follows 

1 



P\ l°° < -He' 



-N 



lL!(iR) 



Li([0,l]) 



< \\u V H1 =E V 

per 



Similarly, 

Using (4.8) we obtain 



< 



d 



eft 



= / p x {t,m) 

J s 

Concerning the term involving arctangents, recalling (1.3) we obtain 



dt 

dt < E u ■ W - s\ 



dt 



2 arctanit^ (t, £(t, x)) 



1 + ui 



uz 



The bounds (4.7) and (4.9) thus yield 

2 arctan^ (s', £(s')) — 2 arctanu^ (s, £(s)) 



< (2||n||^oo + l + 2||P|| L oo) • \s' 

< (lOE u + 1) • \s' - s\ . 



This already provides a bound on the first integral on the right hand side of (4.4). 
Next, call I\,h the last two integrals on the right hand side of (4.4). Notice that 



h+h 



(l + u 2 x {s,y^j -Z v (s';s,y) (l + u 2 x (s', s, y))) 



dy 



(4.6) 



(4.7) 



(4.8) 



(4.9) 



(4.10) 



(4.11) 



Indeed, h+h measures the difference between the measure (l + u 2 (s', y)) dy and the push- forward 
of the measure (l + u 2 (s, x)) dx through the mapping x i— ► s, x). 

Since the push-forward of the measure u 2 dy satisfies the linear conservation law 



w t + (uw) x = , 

comparing (4.12) with (1.5) we deduce 

I u 2 x (s,y) - £ y (s,y)u 2 x (s',£(s';s,y)) dy < [ [ 2\{u 2 - P)u x \dx d 

JO Js JO 

< 2 f (||w||l- + ||P||l-) ||«x||li dt < 2 (4E U + E u ) E u ■ \s' - s\ , 

J s 



(4.12) 



(4.13) 
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because of (4.5), (4.7) and (1.15). Finally, we need to estimate the remaining terms, describing 
by how much the Lebesgue measure fails to be conserved by the transformation x i-> £(s';s,x). 
Observing that 

d -Z y (t,y) = u x (t,Z(t,y))Z v (t,y), ^(0,y) = l, (4.14) 



dt 



we find 



/ I 1 -^y(s';s,y) \dy < / / 

JO Js JO 



[Zv(t,s,y)] 



dy dt 



(4.15) 



: j J £,y( t ; s iy)\ u x{t,£(t;s,y))\dydt . 



To estimate the right hand side of (4.15), we use the decomposition [0, 1] = Y U Y' U Y", where 
Y = {y, Z v (t;s,y) G [(1/2), 2] for alii € [a, s'] } , 

y' = {y, ty(t;s,y) < 1/2 for some t G [s, s'] } , 

Y" = {y; £, y {t; s ,y) > 2 f° r some f g [s, s'] j . 

Integrating over Y one finds 

/ / $ y (t;a,y)\u x (t,Z(t;s,y))\dydt < 2 / ||M*)|| L i < 2£" • |s' - s 

Js JY Js 

Next, if y G V we define 



(4.15) 



r(y) = inf {/ > a; £ y (t;s,y) < 1/2}. 



Observe that y 6 7' implies 



r T (y) 

/ u x {t, (,(t;s,y)) 

J s 



dt > In 2 . 



Therefore 



y/ in 2 Jy/ 



u x (*, £(t;s,y)) 



dt 



dy 



(4.16) 



The estimate for the integral over Y" is entirely analogous, Indeed, the push-forward of the 
Lebesgue measure along characteristic curves from t = s to t = s' satisfies exactly the same 
type of estimates as the pull-back of the Lebesgue measure from t = s' to t = s. All together, 
these three estimates imply 



/ \1-C y (s';s,y)\dy< W E u ■ \s' - s\ 

JYUY'UY" 



(4.17) 



Putting together the estimates (4.6), (4.9), (4.10), (4.3) and (4.17), the distance in (4.4) can 
be estimated by 



(u(s) , u(s')) < 2(1 + E u ) + E u + (10 E* + 1) + 10 {E u f + 10 E* 
This establishes (4.1). 
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s — s . 



(4.18) 
□ 



According to Lemma 5, as long as no peakon interactions occur, the map t >— > u(t) remains 
uniformly Lipschitz continuous w.r.t. our distance functional, with a Lipschitz constant that de- 
pends only on the total energy E u . Since interactions can occur only at isolated times, to obtain 
a global Lipschitz estimate it suffices to show that trajectories are continuous (w.r.t. the distance 
J) also at interaction times. 

Lemma 6. Assume that the multi-peakon solution u(-) contains two or more peakons which interact 
at a time t. Then 



Proof. To fix the ideas, call x = q the place where the interaction occurs, and let pi, . . . ,pk be the 
strengths of the peakons that interact at time r. We here assume that < q < 1. The case where 
two or more groups of peakons interact exactly at the same time r, within the interval [0, 1], can 
be treated similarly. 

For \t — t\ < h, call £ _ (t), £ + (i) respectively the position of the smallest and largest character- 
istic curves passing through the point (r, q), as in (2.5). We observe that u is Lipschitz continuous 
in a neighborhood of each point (r, x), with x ^ q. Hence, for x € [0,1] \ {q} there exists a 
unique characteristic curve t i— » £(t;r,x) passing through x at time r. For a fixed h > 0, the 
transport map ij) is defined as follows. Consider the intervals I_h = [£~(r — h), £ + (r — /i)] and 
Ih = [£~(t + h), £ + (r + h)] . On the complement [0, 1] \ I-h we define 



so that transport is performed along characteristic curves. It now remains to extends ip as a map 
from I-h onto 1^. Toward this goal, we recall that our construction of multi-peakon solutions in 
Section 2 was specifically designed in order to achieve the identity 



lim^ J{u(t - h), u(t + h)) = . 



(4.19) 



i/)(£(T-h;T,x 



■)) =Z(T + h;r,x), 




(4.20) 



For h > we introduce the quantities 





e(h) = 2min {£(-/*), E(h)} 
Notice that (4.20) implies e(h) > and 



max 



{E(-h), E(h)} . 



E{h) -> e (r> g) , 



e(h) -» e (r>$) , 



(4.21) 



as /i — > 0+. Consider the point x* = x*(/i) inside the interval 7-^ = (t — /i), £ + ( r — /i)] , 



implicitly defined by 




I [\ + u 2 x (x)) dx = e(h) . 

C-(t-/») 



For x € [£ (t — h), x*] we define ^(x) as the unique point such that 



/ (1 + u 2 x (t + h, x)) dx = (l + u 2 x (T-h, x))dx. (4.22) 

€-(r+h) h~{r-h) 
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We then extend -0 as an affine map from [x*, £ + (r — h)] onto [tp(x*), £ + (r + h)] , namely 

xl>(e • £+(r - h) + (1 - 6) ■ x*) = 9 ■ £+(r + h) + (1 - 6) ■ 0>(x*) € [0, 1] . 

Finally, we prolong ijj to the whole real line according to (3.3). As usual, the 1-periodic functions 
0i, <p2 are then chosen to be as large as possible, according to (3.4). As h — > 0+, we claim that the 
following quantity approaches zero: 

J^(u(T-h), u(r + h)) 

= [ 1»((x, u(r - h,x), 2, M r - h,x)) , W x). u(r + *, «,)), 2arc t an i ,( T + „, «,)))) 
• 0i(x) (l + u 2 x (t - h, x)) dx 

+ f (l - 0i (x)) (l + u*(t - h, x)) dx 
Jo 

+ [ {l-M^)))(l + ul{r + h^{x)))^{x)dx. 



(4.23) 

It is clear that the restriction of all the above integrals to the complement [0, 1] \ I_h approaches 
zero as h —>■ 0. We now prove that their restriction to I_h also vanishes in the limit. As h — > 0+, 
for i 6 L), we have 

d® ^(x, u(t — h, x), 2 arctan-u^r — h, x)) , (ip(x), u(r + h, ip(x)), 2 arctan u x (r + /t, ip(x)))^J — > 0, 

because all points approach the same limit [q, u(r,q), 7r). The first integral in (4.23) thus ap- 
proaches zero as h — > 0+. 

Concerning the last two integrals, by (4.22) it follows 

01 (x) = 02 (-0(^)) = 1 f° r all z € [£~(r ~ «*] • 

Moreover, our choice of x* implies 

J t {l + u 2 x (T-h, x))dx + J (l + ul(r + h,^(x))j^(x)dx 

< 2max{£(-/i), E{h)} - 2min{£(-h), E{h)} . 

By (4.21), as h — ► 0+ the right hand side of (4.24) approaches zero. Hence the same holds for the 
last two integrals in (4.23). This completes the proof of the lemma. □ 



5 - Continuity w.r.t. the initial data 

We now consider two distinct solutions and study how the distance J(u(t) v(t)) evolves in 
time. To fix the ideas, let t i— > u(t) and t <— ► v(t) be two multi-peakon solutions of (1.1), and assume 
that no interaction occurs within a given time interval [0, T\. In this case, the functions u, v remain 
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Lipschitz continuous. We can thus define the characteristic curves t i-> £(i, y) and t i— > £(i, 2/) as 
the solutions to the Cauchy problems 



(5.1) 

C = t>(*,C), C(0) = y, 

respectively. Let now -0(o) S .F be any transportation plan at time t = 0. For each i G [0, T] we 
can define a transportation plan -0(4) £ F by setting 

1>(t){Z(t,y))±({t, i/> {0) (y)). (5.2) 
The corresponding functions cj)^,^ are then defined according to (3.4), namely 

^(a;) = sup j(9e [0,1]; 6-(l + ul(t,xj) < (l + v 2 x {t, V(t)(x))) i>[ t 



(x) = sup lee [0,1]; 1 + > 0- (l + t£(t, V(t)(^))) ^ (t )(^ 



If initially the point y is mapped into y = tp^(y), then at a later time t > the point £(t,y) 
along the ^-characteristic starting from y is sent to the point £(i, y) along the ^-characteristic start- 
ing from y. We thus transport mass from the point (%(t, y) , u(t,£(t, yj) , 2 arctann^ (t, £(t, y))^j 

to the corresponding point (c(t,y) , f (i, C(*? ^)) : 2 arctanv^ (i, £(i, y))) , 

In the following, our main goal is to provide an upper bound on the time derivative of the 
function 

J^ {t) (u{t) , v(t)) = J d?((x, u(t,x), 2arctanu a; (i,x)) , 

(^(t)(x), v(t, ip{t)(x)), 2 arctan «„.(*, V(t)(z)))) (H-tt^(t,x)) dx 

,i (5.3) 

+ / {l-$\x)){l + u 2 a (t,x))dx 
Jo 



Differentiating the right hand side of (5.3) one obtains several terms, due to 

• Changes in the distance between the points (£, u, 2arctantia;) and (£, v, 2arctant> x ). 

• Changes in the base measures (1 + u^)dx and (1 + v%) dx. 

Throughout the following, by 0(1) we denote a quantity which remains uniformly bounded as u, v 
range in bounded subsets of ifp er . Using the elementary estimate 

\u — v\ < (l + \u\ + \v\) min [\u — v\, 1} , 



19 



we begin by deriving the bound 



I' 1 d 

h = J j t \x-i> (t) (x)\ .$\ x ) (l + u 2 x (t,x)) dx 

< [ u(t,x)-v(t,ip {t) (x)) ■(f>{ t \x)(l + ul{t,x)) dx 
Jo 

<(l+||«(*)|| L - + |Wt)|| L .)-^™(u(t),i;(t)) 
= 0(1)- J*™ (u(t),v(t)). 

Here and in the sequel, the time derivative is computed along characteristics. 
Next, recalling the basic equation (1.1), we consider 



(5.4) 



< 



J ^ u(t,x) - v(t, *p( t )(x)) ■ (f>i\x) (1 + ul(t,x)) dx 



P:(t,x)-PZ{t,i> (t) (x)) .(l + ul(t,x)) dx. 



(5.5) 



In the spatially periodic case, by (1.2) and (2.2) we can write the source terms P™, P% as 



P%(t,x)= l -J X ^{x-y)- u 2 (t,y) + 



u 2 x (t,y) 



dy, 



(5.6) 



■+ / x'(i>( t) (x) - Vw(y)) • vUt ^ t)iy) K ' {t) (y)dy, 

' x— 1 ^ 



where, according to (2.2), 



x'(x) 



e x — e 1 x 
e-1 



< x < 1 , 



X '(x) = X '(x + l) x£M. 



(5.7) 



In the next computation, we use the estimate 

dy = 0(l)-J^(u,v). (5.8) 



J (i + ^(»))-(i + «2(^(y)))^(») 



which holds because of the last two terms in the definition (3.6). Observing that x' is Lipschitz 
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continuous on the open interval ]0, 1[ , we now compute (omitting explicit references to the time t) 
P«(x) - P v x (^(x)) \<\J W(x - y) ■ u\y) - X '^(x) - y) ■ v 2 (y)\ dy 
+ 0(1)- [ x - y - (V»(ar) - ^(y)) 



x-l 



1 

+ 4 



J i X{x - y) (u 2 x {y) - vl(ip(y)) dy 
0(l)-(\x-^(x)\ 



+ \\u 2 - f 2 || L i) 



f 0(1)- I \x-^{x)\- I \y-m\- ^r^ i>'iy)dy 



x-l 



+ 0(l)-^(u,v) + \J^x\x-y)-[^(y)-l] dy 
0(1) ■ \x - i>(x)\ + 0(1) ■ J*(u, v) + 0(1) • (jx - iP(x)\ + £ ^ X "(x - y) ■ ty(y) - y] dy 



= 0(1) • \x - i/>(x) \ + 0(1) • J*(u, v) . 

Integrating over one period we conclude 

h = 0(l)-J^{u{t),v(t)). 
For future use, we observe that a computation entirely similar to (5.9) yields 

P«(x) - P v x {^{x)) | = 0(1) • \x - V(*)| + 0(1) • J*(u,v) . 

Next, we look at the term 



(5.9) 
(5.10) 

(5.11) 



f 1 d 

h = ~t 2arctan u x (t, x) — 2arctan v x (t, V(t)(^)) ■ (f>i\x) (l + u 2 (t,x)) dx . (5.12) 
Jo dt 



Along a characteristic, according to (1.3) one has 

d 2 
— 2arctanu a; (i,£(t)) = — — j 
dt t + u x 

Call 6 U = 2arctanna; , 9 V = 2arctanuj; , so that 

1 o 9 U u T 1 . 



ut 



7T = cos — 

l + ul 2 



1 + ^ 2 



sm( 



— ^^r- = sin 2 — 
l + ul 2 



We now have 



[\l + u x (x))- 
Jo 



u 2 x (x) v x (iP(x)) 



l + ul(x) l + v 2 x (^(x)) 
[\l + u 2 x (x)). ^ 2 ° U{X) - 2 

J 



sin 

< / (l + u 2 (x))- r(x)-^(^(x)) 
Jo 

= 0(1) • J(m,i;) . 
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dx 
v {^{x)) 



dx 



(5.13) 



Next, using (5.11) we compute 

u 2 {x) - P u {x) v 2 (%l>(x)) - P v (i>(x)) 



J 



l + u 2 x (x)) 



l + u 2 (x) 



dx 



< 



l + vl(i>(x)) 

I u 2 (x) — v 2 (il>(x))\dx+ \P u (x) -P v (ip(x)) 
Jo Jo 



dx 



(5.14) 



+ 0(1) 



l + u%(x) l + v 2 {^(x)) 



(l + u 2 x (x))dx 



= 0(1)-J*(u,v), 
where the last term was estimated by observing that 

1 1 



< 2arctanu r — 2arctanu T 



l + u 2 x l + v 2 x 
Putting together all previous estimates we conclude 

Ii + 1 2 + J 3 = 0(1)J*(«,d) 



(5.15) 



(5.16) 



To complete the analysis, we have to consider the terms due to the change in base measures. 
From (1.5) it follows that the production of new mass in the base measures is described by the 
balance laws 

r (i + u% + [u(i + uD] x = [2u 2 + 1 - 2P u \u x = r, 

\ (l + v 2 x ) t +[v(l + v 2 x )] x = [2v 2 + l-2P^]v x = r. 
This leads us to consider two further integrals h,h '■ 

• \2u 2 (x) + 1 - 2P u (x)\ \u x (x)\dx 
= 0(1) ■ J d^ ^(x, u(x), 2 arctan-u^x)) , (i/)(x), v(i/)(x)), 2 arctan 1^(^(2:))^ 
• (l + u 2 (x)) dx 

= 0(1)- J*(u,v). 



< 2 



[2u 2 (x) + l-2P u (x)]u x (x) - 2v 2 (*P(x)) + 1 - 2P v (iP(x)) v x (ip(x)) i(/(x) 
2 {x) -v 2 (i>(x)) + P u (x) -P v (tp(x)) \ \u x (x)\dx 



dx 



(5.19) 



+ 

>o 

V + T" 

J 5 + J 5 



2v 2 (4,{x))+l-2P v (4>(x)) ■ u x (x) -v x (il>(x))il>'(x) 



dx 
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(5.20) 



Using (5.11) we easily obtain 

I's < ^\u 2 {x) -v 2 (^{x))\ + \p u {x) -P"(V(z))|J (l + u 2 x (x))dx 
= 0(1)- J^(u,v). 

On the other hand, recalling (5.13) and using the change of variable y = ip(x), x = i()~ 1 (y), we find 

4' = 0(1)- f u x (x)-v x ^(x))^\x) 
Jo 

u x (x) v x (il)(x)) 



dx 



0(1) 



l + ul(x) l + v 2 {i>(x)) 

+ 0(1)- f v x ^(x))^'(x] 
Jo 



(l + u x (x)) dx 
l + u 2 x (x) 



<0(1)-J^(u,v)+ [ (l + ul(x))-(l + v 2 Mx)))i;'(x) 
Jo 

= 0(1) ■ J(u,v) . 



(5.21) 



dx 



All together, the previous estimates show that 



d 



- J<^> (tt(t), v(t)) <h+h+h + h + I' b + I'l = 0(1) ■ J*n (u(t), v(t)) 



(5.22) 



where 0(1) denotes a quantity which remains uniformly bounded as u, v range on bounded sets of 
Hp eT . As an immediate consequence we obtain 

Lemma 7. Let t u(t), t i-> v(t) be two conservative, spatially periodic multipeakon solutions, 
as in Lemma 2. Then there exists a constant k, depending only on max { ||u||^i er , ||w||i/i er }, such 
that 

J(u{t) , v(t)) < e K| *" s| • J(u{s) , v{s)) s,t€M. (5.23) 



Proof. For t > s the estimate (5.23) follows from (5.22), taking the infimum among all transporta- 
tion plans ip( s ) at time s. The case t < s is obtained simply by observing that the Camassa-Holm 
equations are time-reversible. D 



6 - Proof of the main theorems 

Thanks to the analysis in the previous sections, we now all the ingredients toward a proof of 
Theorems 1 and 2. 

The estimates in (1.16) follow from Lemma 3. Given an initial data u £ Hp eT , to construct the 
solution of the Camassa-Holm equation we consider a sequence of multi-peakons u n , converging 
to u in Hp ei . Then we consider the corresponding solutions t i— > u n (t), defined for all n > 1 and 
t £ JR. This is possible because of Lemmas 1 and 2. 
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We claim that the sequence u n (t) is Cauchy in L 2 ei . Indeed, by Lemma 3 and Lemma 5, 
\u m {t)-u n {t)\\^ < C-J(u m (t), un(t)) < C-e~'*" J{u m (0), u n (0)) < C 2 -e K ^\\u m (t)-u n (t)\\ H1 



Therefore, u n (t) — ► u(t) in L per , for some function u : M i— » H^ er . By interpolation, the conver- 
gence u n — > u also holds in all spaces L£ er , 1 < p < oo. The continuity estimates (1.17)-(1.18) now 
follow by passing to the limit in Lemma 5 and 6. 

It remains to show that the limit function u(-) is actually a solution to the Camassa-Holm 
equation and its energy E(t) in (1.15) is a.e. constant. Toward these goals, we observe that all 
solutions u n are Lipschitz continuous with the same Lipschitz constant, as maps from M into L 2 er . 
Indeed 



*n,a:||L2 er 



+ 



i e -N 

2 



it; 



"7 



(6.1) 



As a consequence, the map t i— » u(t) has uniformly bounded Hp BT norm, and is Lipschitz continuous 
with values in L 2 er . In particular, u is uniformly Holder continuous as a function of t, x and the 
convergence u m (t,x) — > u(t,x) holds uniformly for £ in bounded sets. Moreover, since L 2 is a 



reflexive space, the time derivative u t (t) 6 Lp er is well defined for a.e. t £ M. 



We now observe that, for each n > 1, both sides of the equality 



(6.2) 



are continuous as functions from 1R into L per , and the identity holds at every time t £ M, with 
the exception of the isolated times where a peakon interaction occurs. 

At any time t where no peakon interaction occur in the solution u n , we define to be the 
measure with density u n (t, ■) + \u n ^ x {t, ■) w.r.t. Lebesgue measure. By Lemmas 5 and 4, the map 



(n) 



can be extended by weak continuity to all times t € M. We can now redefine 



P u "(t,x) 



P u (t,x) 



_ e -\x-y\ 



dfit(y)- 



(6.3) 



where fj, t is the weak limit of the measures fi[ n K Because of the convergence J(u n (t), u(t)) 0, 
by Lemma 4 the map t <— > \i t is well defined and continuous w.r.t. the weak topology of measures. 
Using again Lemma 4, we can take the limit of (6.2) as n — > oo, and obtain the identity (1.14), for 
every t G M and P = P u defined by (6.3). 

For each n, the total energy ^ n ' 1 (]0, 1]) = E Un is constant in time and converges to E u as 
n —>■ oo. Therefore we also have 



M*(]0,1]) =i? € 



■/ o 



ti 2 (x) + ti^(x)] 



4 G iR. 



(6.4) 



To complete the proof of Theorem 1, it now only remains to prove that the measure Ht is 
absolutely continuous with density u 2 (t, ■) + ^u^(t, •) w.r.t. Lebesgue measure, for a.e. time t € JR. 



In this direction, we recall that, by (1.5), the function w - 
equation with source 

w t + {uw) x = (u 2 n -P Un )u n>a 



u n satisfies the linear transport 



(6.5) 
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Moreover, along any characteristic curve t by (1.3) one has 



dt 



2 arctan -u n>a; (i,£(i,x)) 



1 + u4 



< 



(6.6) 



whenever x is sufficiently large. For e > small, consider the piecewise affine, 27r-periodic 
function (see fig. 2) 

9 if o < e < 1 

1 if 1 < 9 < tt - e , 

<p{0) = \ (tt - 9)1 e if ^-e<#<^ + e, 

-1 if tt + s<9<2tt-1, 

9-2tt if 2tt-K9 <2tt. 




figure 2 



and define 



0n(t) = / (2 arctan-u^ (t,a;)) u^ x (t,x)dx. 
Jo 

By (1.3) and (6.6) we now have 
^Pn(t) > \- f ul x (t,x)dx-C- [ u^ x dx (6.7) 

al ^ e J {2 arctan u„ ia! e[n— £,7r]U[ — it, — tt+e]} J {2 arctan u n>x 6[ — 1,1]} 

for some constant C, independent of e, n. Since all functions (5 n remain uniformly bounded, by 
(6.7) for any time interval [r, t'\ we obtain 

f f u 2 n:X (t,x)dxdt <eC- (l + r'-r), (6.8) 

J t »/ {2 arctan u„ :a! e[7r — e,7r]U[ — 7r, — 7r+£] } 

where the constant C depends only on the H^ er norm of the functions u n , hence is uniformly valid 
for all n,£. Because of (6.8), the sequence of functions v? nx is equi-integrable on any domain of 
the form [r, r'\ x [0, 1]. Namely 

lim f f u 2 nx (t,x)dxdt = 0, (6.9) 

J{xe[o,i], ' 
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uniformly w.r.t. n. By Lemma 4 we already know that II m£ x (t) — u^.(t) IL j ^ for every fixed 
time t and 1 < p < 2. Thanks to the equi-integrability condition (6.9) we now have 



"7, , — a: in Ij 1 {[t,t'] x [0,1]) 



i 2 - - 2 



By Fubini's theorem, this implies 



lim / u 2 nx (t,x)dx = I ul(t,x)dx 
n ^°°Jo ' Jo 



for a.e. t € [r, r']. At every such time t, the measure /j,t is absolutely continuous and the definition 
(6.3) coincides with (1.2). This completes the proof of Theorem 1. □ 



7 - Uniqueness 

Before proving Theorem 3, we remark that the solution satisfying all conditions in Theorem 
1 need not be unique. 





Example. Let u = u(t, x) be a solution containing exactly two peakons of opposite strengths 
Pi(t) = —p2(t), located at points qi(t) = —q2(t) (see fig. 3). We assume that initially pi(0) > 
0<7i(0) < 0. At a finite time T > 0, the two peakons interact at the origin. In particular, as 
t — > T— there holds 



pi(t) -> oo, p 2 (t) -oo, qx(t) -> 0, q 2 (t) -> . 
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Moreover, ||w(*)||l°° - * 0, while the measure fi t approaches a Dirac mass at the origin. We now 
have various ways to extend the solution beyond time T: 

u(t,x) = 0, (7.1) 

u( T , x) = -u{t - T, -x) , (7.2) 

Clearly, u dissipates all the energy, and does not satisfy the identity (1.15). The function u in (7.2) 
is the one constructed by our algorithm in Section 2. However, there are infinitely many other 
solutions that still satisfy (1.15), for example 

u(t, x) = u(t, x — b) (7-3) 

where u is as in (7.2) and 6/0. The additional condition in Theorem 3 rules out all of them, 
because as r — > T+, the corresponding measures fi T approach a Dirac mass at the point x = b, not 
at the origin. 

We can now give a proof of Theorem 3. As a first step, we extend our distance J to a larger 
domain V, consisting of couples (u,fi), where u £ H^ ev and /j, is a positive (spatially periodic) 
measure whose absolutely continuous part has density u 2 + u 2 x w.r.t. Lebesgue measure. This 
extension is achieved by continuity: 

J((u,n), («,//)) = liminf J(u n ,u n ) (7.4) 

where the infimum is taken over all couple of sequences («„,«„)„>i such that 

\\u n — —> 0, \\u n — — ► , 

U n,x ~^ A* ) ^n,x ~^ A • 

We observe that the flow $ constructed in Theorem 2 can be continuously extended to a locally 
Lipschitz continuous group of transformations on the domain T>. 

Now let t i— > u(t) be a solution of the Cauchy problem (1.1), (1.13), satisfying all the required 
conditions. In particular, the map t i— ► (u(t),fl t ) is Lipschitz continuous w.r.t. the distance J, with 
values in the domain V. 

Calling t i— > (u{t),jit) = ®t(u, u 2 ) the unique solution of the Cauchy problem obtained as 
limit of multi-peakon approximations, we need to show that u(t) = u(t) for all t. To fix the ideas, 
let t > 0. By the Lipschitz continuity of the flow, we can use the error estimate 

j((fi(t),£t), (u(t),A*t)) < e C2t j liminf J ({u(t + h) , jl T+h ) , $ h (u(r), jl T )^ dr (7.5) 

J 

For a proof of (7.5), see pp. 25-27 in [B]. The conditions stated in Theorem 1 now imply that, at 
almost every time r, the measure jl t is absolutely continuous and the integrand in (7.5) vanishes. 
Therefore u(t) = u(t) for all t. □ 
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We can now prove that, in multi-peakon solutions, interactions involving exactly two peakons 
are the only possible ones. 



Corollary. Let t u(t,) be a multi-peakon solution of the form (1.7), which remains regular on 
the open interval ]0, T[ . Assume that at time T > an interaction occurs, say among the first k 
peakons, so that 

lim qi(t) = q i = l,...,k. 

Then k = 2. 

Proof. We first observe that the Camassa-Holm equations (1.1) are time-reversible. In particular, 
our proof of Theorem 3 shows that the solution to a Cauchy problem is unique both forward and 
backward in time. 

Now consider the data (u(T), ht) £ where \xt is the weak limit of the measures Ht having 
density u 2 (t) + u 2 (t) w.r.t. Lebesgue measure, as t — > T-. By the analysis in Section 2, we can 
construct a backward solution of this Cauchy problem in terms of exactly two incoming peakons. 
By uniqueness, this must coincide with the given solution u(-) for all t € [0, T\. □ 
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